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Abstract 



Global aspects of the motion of passive scalars in time-dependent 
incompressible fluid flows are well described by volume-preserving (Li- 
ouvillian) three-dimensional maps. In this paper the possible invari- 
ant structures in Liouvillian maps and the two most interesting nearly- 
integrable cases are investigated. In addition, the fundamental role of 
invariant lines in organizing the dynamics of this type of system is ex- 
posed. Bifurcations involving the destruction of some invariant lines and 
tubes and the creation of new ones are described in detail. 
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1. Introduction 

The dynamics of deterministic passive scalars suspended in incompressible fluid 
flows has very different characteristics as the dimensionality of the flow varies. In 
the most general formulation of the problem, we are interested the study of the 
trajectories given by the following set of differential equations: 

r = u(a;,y,z,t), (1) 

where the velocity field u(a;, j/, 2, satisfies 

V • u = 0. (2) 

Several simpler instances of the problem have been investigated thoroughly in the 
past decade/ Let us quickly review some of them. 

When the flow is two-dimensional and stationary, Eq. (1) implies that the 
velocity field u(a;, y) can be derived, from a stream function ^'(a;, y) 

«. = ^;% = -^. (3) 

In this case the equations of motion have a Hamiltonian structure with \1/ playing 
the role of the Hamiltonian. Hamiltonian systems with only one degree of freedom 
are always integrable. For bounded systems, this means that all the trajectories in 
the phase space are closed curves. Accordingly, trajectories of passive scalars in 2D 
steady flows will have the same topology. 

Next in order of increasing complexity comes the class of time-periodic 2D 
flows. ^ Here the stream function becomes time dependent and the system is equiv- 
alent to a non-autonomous Hamiltonian one. In this case, only integrable systems 
show a phase space completely foliated into invariant tori as in the previous case. 
However, the KAM theorem states that a measurable portion of these invariant tori 
persist in phase space even for flnite departures from integrability. These invariant 
surfaces are relevant for passive scalar applications since they represent dynamical 
barriers for transport across space. 
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Steady 3D flows are somewhat equivalent to the case described in the previous 
paragraph. Arnold has shown that the motion of passive scalars in steady inviscid 
fluid flows is integrable if = V x u is almost nowhere parallel to u. On the 
contrary, if the Beltrami condition a; = Au is satisfied, it has been found that 
trajectories which are chaotic coexist with regular ones. The regular trajectories 
cover two-dimensional surfaces which are invariant under the time evolution of the 
passive scalar. In turn, the invariant surfaces separate the available space into 
disconnected regions which confine the chaotic trajectories. 

Next in complexity is the case of 3D flows with periodic time dependence. In 
this case, the form of the equations of motion is no longer Hamiltonian. However, the 
stroboscopic maps of such systems are three-dimensional and, by virtue of Eq. (2), 
volume preserving. In the following we will refer to these maps as 3D-Liouvillian 
or in short, 3DLM. The invariant structures displayed by these maps in the nearly 
integrable regime have been investigated in detail in the past^ because of their 
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relevance to the mixing efficiency ' of fluid flows. In particular, we found one class 
of 3DLM where these structures divide the space into separate regions, preventing 
a single deterministic passive scalar trajectory from reaching the whole space. In 
another class of 3DLM, however, unbounded deterministic diffusion has been found 
arbitrarily near the corresponding integrable map. 

We have also conjectured that the objects organizing the dynamics of 3DLMs 
are one dimensional quasiperiodic orbits which we name invariant lines. In this 
paper we review the former results and provide further evidence to support this 
picture. 

Sec. 2 is devoted to the nearly-integrable behaviour of the 3DLM referred to 
above. In Sec. 3 we study the evolution of the structures related to invariant lines 
and tubes around them as the parameters of the maps are varied. We show that the 
disintegration of an invariant line into flxed points or periodic orbits at parameters 
values that make the line strongly resonant is accompanied by the birth of a new 
invariant line. The ubiquity of invariant lines is shown to follow from the generic 
topology of the invariant manifolds of the periodic orbits. Finally, in the concluding 
Sec. 4 we summarize our results and discuss some problems for future research. 
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2. Three-dimensional maps and their invariant structures 

We concentrate our attention on volume-preserving maps L : — > on the 
three-dimensional torus of the form 

{x', y', z') = {x^ u{x, z),y + v{x', z), z + w{x', y')), (4) 

where u, v, and w are doubly-periodic functions. In particular, we shall investigate 
the following truncation of the Fourier expansions for these functions: 

x' = X + Ai sin z + C2 cos y (mod 2tv), 
y' = y + Bi sin x' + A2 cos z (mod 27r), 

z' = z + CJ{y') + B^gix') (mod27r), (5) 

where the functions / and g are smooth on the circle. 

In general we call integrable a map L/,^ for which a set of variables I £ R*', 
9 e T^~'^ can be found such that 

1=1, 

e' = e + uj{i). (6) 

Clearly, the variables I, which we call actions by analogy with Hamiltonian dy- 
namics, parametrize a family of (3 — /c)-dimensional geometrical objects that are 
invariant under the dynamics. In turn, the 9 variables specify the position of a 
point on a given invariant object. At each iteration the angular variables rotate by 
an angle given by the corresponding component of the (3 — /c)-dimensional vector 
uj. This angle might in general vary from one invariant set to another, but remains 
constant for all the points on each of these sets. Let us describe the four possible 
cases: 

(a) k = 0. Since there are no actions, the frequency vector co is constant and 
this integrable case corresponds to a uniform rotation on T^. 

(b) k = 1. The motion takes place on two-dimensional tori defined by / = 
const, and a point on one of these tori is specified by the values of two angles 
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9i,92. Therefore, the motion on these two-tori is a uniform rotation for which the 
frequency depends only on the value of the action. 

(c) k = 2. In this case, there are two actions /i, I2 which parametrize a family 
of invariant circles (one-dimensional tori). The angular variable rotates with a 
frequency c<;(/i,/2). 

(d) k = 3. In this case, the map in Eq. (6) is the identity and each point of 
the space (parametrized by /i,/2 and I3) is invariant. There are no angles in this 
problem. 

For case (a) the results of Ruelle and Takens^'^'^'^ indicate that small perturba- 
tions of the integrable system will produce completely chaotic maps. Although the 
implication of volume preservation for these results is a very interesting problem, it 
lies outside the scope of this work, which aims towards the understanding of lower- 
dimensional invariant objects. At the opposite extreme, /c = 3, the nearly integrable 
systems are small perturbations of the identity I' = I -|- eF(I) or AI = eF(I). In 
the e ^ limit, the dynamics of this system can be described by a set of three 
autonomous differential equations after an appropriate rescaling of the time. This 
case (d) also lies outside the interest of the present paper. Therefore, we shall con- 
sider only small perturbations Li and L2 around the respective integrable maps 
and L2. 

An example of nearly-integrable maps with a single action is the map Li given 
by the following restriction of Eq. (5): 

x' = X + A sin z + eac cos y 

y' = y + eas sin x + A cos z 

z' = z + eac siny' + eaB cos x' . (7) 

At e = this map is integrable with only one action variable z. The motion in this 
limit takes place on the surfaces of constant z, and the angular variables x and y ro- 
tate at each iteration by constant angles lVx = A sin z and ujy = A cos z respectively. 
The map in Eq. (7) corresponds to a discretization of the ABC model proposed by 
Arnold^ as an example of chaotic streamlines in a stationary 3-dimensional flow, 
and extensively studied by Dombre et al. The discretization is equivalent to the 
addition of periodic time dependence to the flow. 
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In order to investigate the effect of a small nonintegrable perturbation of am- 
plitude e on the form of the invariant surfaces, we implement a perturbative scheme. 
The condition for a surface defined by the equation r = r(t, s) (where r = (a;, y, z) 
and t, s are parameters) to be invariant is 

U{v{t,s)) = v{t\s'). (8) 

We start with the invariant surface 2; = 2:0 of the integrable case and suppose that 
the perturbed surface can be written in the form z = zq -\- '^e^Hn{x,y). Here 
the parameters t and s are identified with the coordinates x and y. Inserting 
this expression into the invariance condition, we obtain an infinite system of linear 
functional equations for the unknown functions that in principle can be solved 
order by order. The first step in this expansion leads to 

z' = Zq + eHi{x' , y') + O(e^) — zq -\- eHi(x + Asinzo, y + Acoszq) + O(e^) 
= zo + eHi{x,y) + eac sin(y + Acoszq) + ecus cos(a; + Asinzo) + 0{e^). (9) 

Consequently, Hi{x,y) satisfies the linear functional equation 

Hi{x + Asi'nzo, y + Acoszo) = Hi{x, y) +ac sin(^/ + A cos zoj+as cos(a; + Asinzo). 

(10) 

The previous equation can be easily solved by expanding Hi (a;, y) in a double Fourier 
series. Setting Hi{x,y) — Yli '^mn&^^^^'^'^^^ in (10) we obtain 

r i(mAsinzo+nAcos2o) 1 1 _ \X , p-iAcos zo _ s: McoszolX 

I IX ^iAsinzo,: —iAsinzols: 

+ -^[t>l,me + ()-l,me \<JO,n- [i-i-) 

Notice that the coefficients amn in Eq- (H) remain undefined if the resonance con- 
dition 

mA sin zq + nA cos zq = 2'Kk (12) 

is satisfied for {m,n,k) — (±1, 0, A;), (0, ±1, /c). However, they can be uniquely 
determined by requiring that amnS be continuous functions of z. Thus, the only four 
non-zero coefficients in Hi{x,y) are ao,±i and a±i,o- This is, however, a peculiarity 



-7- 



emerging from the form of the nonhnear terms in the equation for the action z. If 
higher-order Fourier components were considered in the original map, then more 
non-zero coefficients would appear in the expression for H\. Also, in our case, when 
the expansion is carried over to higher orders in e, the corrections H^^ will contain 
correspondingly-highcr-order Fourier coefficients. 

Eq. (12) has stronger consequences when it holds for {m,n,k) = (±1,0, /c) 
and {m,n,k) = (0, ±l,/c). Namely, the values of ao,±i and a±i^o respectively are 
diverging. This occurs, of course, only for few special values of zq. The invariant 
surfaces where the condition (12) is satisfied, are called resonant. It is clear from 
the perturbative arguments that these resonant surfaces will exhibit a strongly- 
singular behaviour at finite e. To understand this special behaviour notice first that 
integrable motion on a resonant surface is such that each individual trajectory is 
not dense on the surface but rather fills an invariant curve contained in it. For e = 
there is a continuous family of such invariant lines covering the entire surface. At 
e 7^ however, only a finite (and even) number of such lines survives the presence 
of the nonlinear perturbation. This occurs by a mechanism similar to the Poincare- 
Birkhoff phenomenon for periodic orbits in 2D. Moreover, as in the latter case, it 
turns out that half of these lines are stable (elliptic) and the other half, unstable 
(hyperbolic). In the next section we will illustrate the origin of such behaviour 
with a perturbative calculation. At non-zero perturbation strength e a family of 
elliptic invariant tubes is formed around the stable lines. Similarly, associated with 
the unstable lines an H-shaped chaotic slab emerges. This scheme is repeated for 
higher-order resonances at the corresponding order of the perturbation expansion. 

For the surfaces where the frequencies u are far from satisfying the resonance 
condition, we might expect that some version of the KAM theorem will hold and 
that slightly-deformed invariant surfaces will persist if e is not too big. Actually, 
there is extensive numerical evidence in favour of such a conjecture. In Fig. 1, the 
various types of trajectories are shown for the map of Eq. (7). 

By comparing this picture with the analogous one for 2D-Hamiltonian systems, 
one is tempted to conjecture that the organizing role played by the periodic orbits 
in the latter should be assigned here to the invariant lines. In fact, in the same way 
that invariant tori of the 2D case can be systematically approached with sequences 
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of periodic orbits, in our case a similar (but unfortunately not so systematic) ap- 
proximating strategy can be designed by using invariant lines instead. In Sec. 3, 
we will see that one-dimensional invariant objects pervade the phase space of these 
systems and their presence seems to be a robust property. 

The reader can easily recognize that the preserved invariant surfaces are barri- 
ers through which a chaotically moving particle cannot penetrate. In other words, 
diffusion of an individual trajectory throughout the entire space is not allowed. We 
will show, however, that the opposite holds for L2 maps. 

Let us now consider the maps close to the integrable case with two almost- 
conserved quantities, L2: 



At e = the motion takes place on lines I = (Ii, I2) = const. We could examine the 
perturbative behaviour of those invariant lines with the hope of finding some sort 
of KAM result for them. However, by a mechanism similar to the one responsible 
for breaking the resonant surfaces in Li maps, all the invariant lines are destroyed 
to first order in e. Later, we shall illustrate this process with a particular example. 
In fact, we can understand the origin of the singularity of this integrable case by 
means of the following argument based on an adiabatic approximation. Suppose 
that c<;(I) is irrational for some given values of the arguments. In the limit of e — 0, 
we can assume that before I changes significantly, the angle 9 covers uniformly the 
entire (0, 2it) interval. Under these circumstances, the variation of I can only be 
sensitive to averages of P = (Pi, P2) over all the possible values of 9. Therefore 



where < >$ stands for the ^-average. Thus, the dynamics of the action variables 
decouples from 9 for non-resonant oj. Eq. (14) leads in the limit e ^ to a system 
of two ordinary differential equations 



I[ = Ii + ePiih, 9)- I!,=l2 + eP2{I[, 9) 



9' = 9 + u;iI[J^). 



(13) 



AI = e<P(I,^) >0=eP(I) 



(14) 



P(I) 



(15) 
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where the identification e = At was made. This system can be easily integrated by 
variable separation. Its trajectories l{t) satisfy 

f ' Pi{I)dI - [ ' P2{I)dI = W{h,l2) = const. (16) 
Jo Jo 

In other words, in the e — > limit, the action variables slowly evolve along the 
curves defined by Eq. (16). Including the fast motion in the 9 direction, we infer 
that the originally-invariant lines parallel to the 6 axis coalesce in invariant surfaces 
"Ep defined by the condition W{Ii,l2) = p. A typical trajectory will densely cover 
such surfaces rather than move on an invariant curve. The adiabatic approximation 
is exact in the limit e — 0. One is then tempted to conclude that the situation 
is similar to the one described for one-action maps in the sense that the adiabatic 
invariant surfaces can survive when finite nonlinearitics arc present. However, a 
new and very interesting phenomenon appears in L2 maps. We first notice that 
the adiabatic approximation is bound to fail whenever the resonance condition 
a;(/i, I2) = 2nk/n is satisfied. Moreover, this condition defines a family of surfaces 
which in a generic case does not coincide with the family of invariant surfaces. As 
a consequence, each invariant surface will intersect at least one resonance sheet. At 
the intersections, the adiabatic approximation is spoiled and so is the smoothness 
of the invariant surfaces. Far from the intersection, however, one expects that 
trajectories evolve on surfaces which are slight deformations of the ones given by 
Eq. (16). In order to exemplify the characteristic behaviour of L2 maps, we use an 
appropriate restriction of the family of maps defined in Eq. (5): 

x' = X + eaAi sin z + eaca cos y 
y' = y + eaBi sin x' + eaA2 cos z 

z' = z + C^f{y') + B2g{x'). (17) 

For e = 0, the map in Eq (17) is integrable and has the form of Eq. (13). In this 
case the lines corresponding to constant values of x and y are invariant. A search 
for invariant lines of the form (x, y) = {xq, yo) + '^{Xn{z), Yn{z))€^, in the nearly- 
integrable case can be performed perturbatively. The order-e calculation leads to 
the following couple of functional equations 

X[z + Cif{yo) + B2g{xo)] = X{z) -I- a^i sin2; -|- cosyo, (18a) 
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Y[z + Cif{yo) + B2g{xo)] = Y{z) + as^ sin xq + qjaj cos z. 



(186) 



Clearly, the constant terms containing xq and j/o in Eq. (18) lead to the diver- 
gence of the zeroth-order Fourier coefficient of the X{z) and Y{z) functions. Since 
< sin 2; >z=< cos 2; >z= 0, the adiabatic invariant surfaces of Eq. (16) become 



Fig. 2 shows some of these surfaces projected down to the z = plane. To 
illustrate the effect of the resonances, we plot in Fig. 3 a trajectory of the map in 
Eq. (17) for a non-zero but very small value of e. The location of the lowest-order 
resonance a; = for a particular election of the functions / and g is indicated by 
the dashed line. Close to this line, the trajectory oscillates wildly and jumps from 
one adiabatic surface to another. 

To depict this behaviour in a different way, we show in Fig. 4 the time evolution 
of the value of Wq, which would be constant if the adiabatic approximation were 
exact. One can easily recognize intervals where Wq is almost constant followed 
by relatively short periods of oscillatory behaviour. Naturally, these oscillations 
occur whenever the trajectory crosses the first-order resonance. Notice that as a 
consequence of these oscillations, Wq randomly jumps from one asymptotically- 
constant value to another, corresponding to two different adiabatic surfaces. 

One striking consequence of this dynamical behaviour is that a single trajec- 
tory can in principle visit the entire region of space where the adiabatic surfaces 
which intersect with the first-order resonance reside. The size of this region can be 
controlled by choosing the functional form of the frequency oj. In Fig. 5. we show 
two nearly-extreme cases. In the first one (Fig. 5a), a case where the resonance 
condition almost coincides with one of the adiabatic invariant surfaces is shown. 
Several initial conditions have been used to generate this picture. Notice that most 
of the trajectories evolve on smooth surfaces which are roughly the same as those 
shown in Fig. 2. In addition, one can observe a small region of chaotic trajec- 
tories associated with the invariant surfaces intersecting the first-order resonance 
condition. Of course, higher-order resonances have similar effects in other regions 
of the space, but these are not evident on the time scale of the picture. On the 



Wo{x,y) — sinj/ -|- ckbi cos a; = /3. 



(19) 
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other hand, Fig. 5b shows the opposite extreme. Here, the first-order resonance 
indicated by the dashed fine intersects almost all the surfaces of Eq. (19). In this 
figure, the iterations of only one trajectory are shown. It is now apparent that 
this single trajectory visits all the available space. The rate of diffusion, D, has 
been estimated analytically in Ref. 12. For the case of Fig. 5b it is shown that 
D = 0{e^). Numerically, we found that D = 0{e"') where 7 = 2.0 ± 0.3. 

3. Tubes and invariant lines 

The standard approach to the systematic study of maps and their properties is 
based on the principle of decomposition in terms of the simplest type of invariant 
objects. Since in 3DLM, periodic orbits arc gcnerically unstable, the lowest- 
order invariant objects which underlie the organization of the dynamics are the 
invariant lines. In the case of Li maps, invariant lines behave in a way which is 
reminiscent of the fixed points in 2D conservative maps. In particular, slices through 
the dynamics on the (x, y, 2)-torus transversal to one of the angle directions, x or 
y, are similar to phase-space portraits of the standard map for example. To some 
extent however, this similarity is misleading. It gives the wrong impression that 
a simple extrapolation to one more dimensions of the theory for periodic orbits is 
sufficient for the understanding of invariant lines in 3DLM. In fact, while periodic 
orbits are solutions of algebraic equations, in order to find invariant lines one has 
to solve functional equations. In this section, we shall illustrate this difference 
through the study of the simplest bifurcations which the invariant lines undergo. 
By analogy with the case of periodic orbits, bifurcations are qualitative changes in 
the properties of invariant lines occurring as the parameters of the map vary. In 
what follows, wc restrict ourselves to bifurcations in which invariant lines are either 
created or destroyed. 

The simplest and most dramatic destruction of invariant lines occurs when the 
nonintegrability parameter, e, is turned on. As described in Sec. 2, the continuous 
families of invariant lines present in the integrable case are replaced with a few stable 
tubes and a chaotic slab (see Fig. 1). To illustrate this process we can implement 
a perturbative scheme for invariant lines which is similar to the one already used 
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for surfaces. For the sake of clarity we restrict ourselves to the (±1, 0, 0) resonant 
surfaces. In the integrable case this resonance occurs if Ux = ^sin2;o = i.e., at 
zq = Zi with zi = and Z2 — tt. It is clear that the trajectories on these resonant 
surfaces lie on lines parallel to the y-axis. Each initial condition xq on the resonant 
surface corresponds to a different line which is consequently defined by a; = a;o 
and z = Zi- We want to understand what happens with these lines when e ^ 0. 
Therefore, we will look for perturbed lines of the form 

x = xo + eX(y) + 0{e^) ; z = Zi + eZ{y) + Oie^) (20) 

and confine our calculations to order e. By requiring invariancc of (20) under the 
iteration of the map, we find two functional equations, for X{y) and Z{y) 

X{y + A cos Zi) = X{y) + AZ{y) cos Zi + cosy, (21a) 
Z(y + Acoszi) = Z{y) + sin(y + A cos Zi) + cosxq. (21^*) 

As before, expanding the two functions in Fourier series, X{y) — ^a„^e*"^^ and 
Z{y) — ^a^^e*"^ we obtain for the coefficients 

^^x^ginAcosz, Acoszian' + ^(5l,n + 5-1,^), (22a) 

^^.(ginAcos., ^[^_^^^g-Mcos., _ g^^^^iAcos z,^ ^ 5o,n«B COS O^q. (226) 

Notice that when n = 0, the factor (e^'^^cos^i _ vanishes whereas the r.h.s of 
Eq. (22b) does not. Therefore, the first order correction for the invariant lines is 
finite only if cosxq = 0. This implies that out of the infinity of invariant lines 
corresponding to the (ihl, 0, 0) resonant surfaces zq = Zi, only those defined by the 
(xo,Zi) pairs (7r/2,7r), (37r/2,0), (7r/2,0) and (37r/2,7r) survive at first order in e. 
Numerical computations indicate that half of these lines are dynamically stable and 
the other half are unstable. This result is a remarkable manifestation in 3DLM of 
a scenario similar to the one which, for Hamiltonian systems, is predicted by the 
Poincare-Birkhoff theorem. While around stable lines, a family of elliptic invariant 
tubes is formed, associated with the unstable lines an H-shaped chaotic slab emerges. 
This scheme is repeated for higher-order resonances at the corresponding order 
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of the perturbation expansion. The location of the lowest-order elliptic lines is 
indicated schematically in Fig. 6. 

A different bifurcation, in which tubes are created rather then destroyed, takes 
place dX A = 2n. For A > 27r, one has additional 0(e) solutions to the invariant- 
surface resonance condition of Eq. (12) of the type (m,n, /c) = (±1,0, ±1) and 
(m, 71, fc) = (0, ±1,±1). As before, these solutions correspond to elliptic tubes 
embedded in a chaotic slab. Since all the new tubes have similar properties, we 
only consider the (0, ±1,±1) ones. For these tubes, Ux = ±V A'^ — 4%'^, u)y = 2ti 
for all values of A, and therefore, for A > 27r, they are parallel to the x-axis. 
Their average position in z, Zi, satisfies cos^^ = ±27rA~^. Unlike the tubes of 
Fig. 6, the new tubes travel along the z-direction as the parameter A changes. 
Accordingly, we shall refer to them as travelling tubes. At A = 27r, the travelling 
tubes are created at Zi = 0, tt, at apparently the same position at which the old 
(±1, 0, 0)-tubes lie in the perpendicular direction. One would naively expect that 
the collision between the new and the old tubes would lead to an increase in the 
degree of chaos in this part of space. While numerical investigations^'^ do indeed 
confirm this expectation, it turns out that the actual scenario is more subtle. In 
fact, at A = 2n, the (±1,0,0) invariant lines themselves become resonant (see Eq. 
(22)). This is simply a manifestation of the fact that at e = 0, each of these lines 
is a continuous family of fixed points with Uy = 27r. Accordingly, at finite e, the 
corresponding tubes degenerate into chaotic trajectories in a range of 0{^/e) around 
A — 2n. Remarkably, it is within this range that the travelling tubes are born. 
These appear as regular trajectories which glue together pairs of fixed points of the 
type {uJx,ujy,uJz) = (0,27r, 0). Such fixed points exist in the interval A G {2tv — 
eaBdc-i 2'K+eaBCic) at (xq, 2/o, ^o) such that, to lowest order in e and with 5 = A—2tt 



sin2;o = ± — \/e'^aQ{l — a^) + 5'^. 

27T * 

The ± signs in Eq. (23) are restricted such that only eight fixed points are obtained 
which, at the lower end of the existence interval, are pairwise degenerate. Here, 




(23) 
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the pairs lie at (f , 0, tt), (f , tt, tt), 0, 0) and tt, 0). For A > 27r - easac 
this degeneracy is lifted and the fixed points belonging to a pair drift way from 
each other mainly in the a;-direction. For example, at ^ = 27r the members of each 
pair are two quadrants apart in x. Moreover, at the upper end of the existence 
interval, the fixed points in each pair collide after having wrapped once around the 
a;-direction. 

We now turn to discussing the properties of trajectories that lie in the vicinity 
of the fixed points. For this, one needs to understand the behaviour of the corre- 
sponding stable and unstable manifolds. The fixed points in each pair have two 
complex-conjugate eigenvalues and one real. While, for one of the fixed points, the 
complex eigenvalues correspond to a 2D-stable manifold and the real eigenvalue 
to a ID-unstable manifold, for the second fixed point it is the other way around. 
Moreover, when A is just slightly larger than 27r — easac, the manifolds correspond- 
ing to the real eigenvalue are directed towards the second fixed point of the pair. 
However, gcncrically the two ID manifolds do not meet but rather wind up around 
each other approximately filling up a tube which runs between the two fixed points 
and then opens up along the complex manifolds. This structure of the manifolds 
generates a region of chaotic motion that can be considered as the generalization 
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of Shilnikov chaos to 3D maps. Remarkably, this scenario parallels the one in 2D 
Hamiltonian maps. Notice however that while in 2D the intersections of the man- 
ifolds are unavoidable, in 3D such intersections are non-generic. The parallelism 
between the two can instructively be pursued further. In the same way as in 2D, 
the heteroclinic chaotic regions are usually bounded by regular quasiperiodic mo- 
tion on invariant curves encircling the elliptic periodic orbits; in our 3D example, 
the Shilnikov-type chaos is bounded by a nested family of toroidal surfaces, each 
hosting three- frequency quasiperiodic motion. Through the interior of this family 
runs a circular invariant line (which comes as a replacement for the elliptic point in 
the 2D picture) and it is threaded by both the stable/unstable manifold pair of the 
fixed points and their associated Shilnikov-chaotic trajectories. A trajectory lying 
on one of these doughnut-shaped surfaces is shown in Fig. 7. 

As the fixed points move apart, they pull the tori along, stretching them into 
travelling tubes by the time they collide on the other side of the a;-circle (see Fig. 



-15- 



8). One expects similar bifurcations leading to new travelling tubes to take place 
for all A = 27rZ with Z = 1, 2, Therefore, at large enough A, space will be mostly- 
filled with such tubes rather than with invariant surfaces as is the case for A <2'k. 



4. Concluding remarks 

Three-dimensional Liouvillian maps present an extremely rich variety of dy- 
namical phenomena. Many of these phenomena show an analogy in the nearly- 
integrable limit to the behaviour of Hamiltonian systems of either two or three 
degrees of freedom, which generate symplectic maps in two and four dimensions 
respectively. The similarity with one or the other is governed by the number of 
invariant quantities of the corresponding integrable system. 

Recent progress has been made in the direction of extending KAM results to 
this kind of dynamical system. However, there are several interesting problems 
which are still open and whose solutions, although probably relying on straight- 
forward extensions of the approaches used in other dynamical systems, require the 
development of non-trivial techniques. For example, in the case of maps with one 
action, the smooth two-dimensional KAM surfaces break at a critical value of the 
nonlinearity parameter. The question of how these surfaces behave at the break- 
down point arises as a natural extension of similar studies in two-dimensional area- 
preserving maps. In our context, however, the problem becomes a version of the 
still-unsolved transition to chaos in three-frequency systems. Our studies strongly 
suggest that its solution requires the ability to investigate the invariant lines playing 
the role of periodic orbits in the lower-dimensional case. In fact, we have shown in 
Sec. 3 that the appearance of these objects is generic even in the parameter regions 
where the, in this case more elusive, periodic orbits are present. 

The ergodic properties of L2 maps should be extensively investigated. We were 
able to obtain good estimates for the diffusion rate in the regions dominated by the 
first-order resonances. However, one can imagine a situation where a fraction of the 
adiabatic surfaces do not cross any first-order resonance. In these regions the local 
diffusion rate is determined by the resonances with n > 2. A better understanding 
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ol the interaction between the adiabatic and resonant motions is required to obtain 
a similarly-rehable estimate for this case. 

Finally, it is important to remark that the codimension of the whole class of 
Liouvillian-maps is higher than that of the families that we have studied. Therefore, 
it is very relevant to applications in hydrodynamics to investigate (a) the conditions 
for which real 3D time-periodic flows approach the subclass presented here, and (b) 
to what extent the reported behaviour extends beyond the limits of such a subclass. 
Research in these directions is currently in progress. 
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Figure Captions 

Fig. 1. Numerically obtained invariant surfaces and chaotic volumes. The parameter 
values in Eq. (7) are A = 1.5, as = 1, ac = 2 and e = 0.1. (a) In order 
of increasing z: non-resonant surface {xq = 0, yo = 0.56, zq = 0.11) and 
tubular surfaces around the (1, —1, 0), (1, 0, 0) and (0, 1, 0) resonances -{xq = 0, 
yo = 0.9, zo = 0.395), {xq = 0.75, yo ^ 0, zq ^ 0.4626) and {xq = 0, yo = 0.5, 
zq — 0.828) respectively, (b) H-shaped chaotic volume associated with the 
hyperbolic line of the (0, 1, 0) resonance {xq — 0, yo — 0.5 and zq — 0.24). All 
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the initial conditions are indicated in fractions of 27r and the box represents 
the [0, 2txY region. 

Fig. 2. A few members of the family of surfaces defined by Eq. (17) for ctSi = 1-5 and 
= 2 projected down to the z = Q plane. For comparison, the dashed curve 

is the first turn of the trajectory depicted in Fig. 3. 
Fig. 3. One trajectory of the map in Eq. (17) for = 1, — 2-5, Ci = i?2 = 4 

and e = 0.001, / = cos and g = sin. and are the same as in Fig. 2. 

The dashed line indicates the location of the lowest-order resonance. 

Fig. 4. Wo{xn,yn) (Eq (19)) vs. n. The map and parameter values are the same as in 
Fig. 3. 

Fig. 5. Iterations of the map in Eq. (17) lying in the slice < z < 0.01. The dashed 
lines indicate the location of the lowest-order resonances, (a) Ci = 2.5, B2 = 4, 
/ = sin, g = cos and the remaining parameters are as in Fig. 4. Several initial 
conditions distributed along the j/-axis and the x = tt line were necessary to 
obtain this picture, (b) /, g and all the parameters are the same as in Fig. 3. 
Only one initial condition is required to generate this picture. 

Fig. 6. Schematic representation of the lowest-order elliptic resonances. The dashed 
lines parallel to the x and y-axis represent the stable invariant lines of the 
(0,1,0) and (1,0,0) resonant surfaces respectively, which survive when e 7^ 
in Eq. (7). The cylinders sketch the tubular invariant surfaces that appear 
around these lines at finite e. 

Fig. 7. A trajectory of (7) which lies on a toroidal invariant surface. Here, A = 
6.276, CKs = l,Q!c = 2 and e = 0.01. The pluses denote the pair of fixed 
points associated with this trajectory. The coordinates are expressed in units 
of 27r. 

Fig. 8. A brand new travelling tube . The parameters are the same as in Fig. 7 only 
that ^ = 6.3. Here as well as in Fig. 7, the coordinates are expressed in units 
of 27r. 



